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Abstract. We present an algorithmic way of exactly computing Belyi functions for 
hypermaps and triangulations in genus 0 or 1, and the associated dessins, based on a 
numerical iterative approach initialized from a circle packing combined with subsequent 
lattice reduction. The main advantage compared to previous methods is that it is applicable 
to much larger graphs; we use very little algebraic geometry, and aim for this paper to be as 
self-contained as possible. 
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Introduction 

Let T be a triangulation of the sphere § 2 , i.e. a planar graph embedded in the sphere in 
such a way that all the faces it delimits have degree 3. From T, one can construct a complex 
structure on the sphere by gluing together equilateral triangles according to the combinatorics 
of T ; this leads to a manifold with conical singularities at the locations of the vertices of T, 
and the Riemann uniformization theorem ensures the existence of a conformal bijection dq- 
between this surface and the Riemann sphere C*. The image of T by <lv is a triangulation 
of the Riemann sphere with the same combinatorics as T and edges embedded as analytic 
curves, which is well-defined up to the action of the Mobius group. 

While this construction is quite simple, computing the embedding explicitly is very difficult 
in general, as soon as the number of vertices of T is not very small. One of the goals of this 
paper is to describe an algorithm to do it automatically, in a “semi-numerical” way (that is, 
going through an approximate, numerical stage to obtain an exact, explicit outcome). 

Assume in addition that all the vertices of T have even degree. This is equivalent to two 
other properties of the triangulation: first, its faces can be partitioned into two sets in such 
a way that a face in one set is only adjacent to faces in the other set; this is usually done 
graphically by coloring one half of the faces white and the other half black. Second, the 
vertices of T can be partitioned into three classes in such a way that along the boundary of 
each face, exactly one vertex of each class appears. 

Choosing the triangles used to build our Riemann surface to be the two hemispheres of 
the Riemann sphere separated by the extended real line, with vertices at 0, 1 and oo, from 
the data of T we obtain a covering ttt : C* -A C* that is ramified only above {0,1, oo} (such 
a map is known as a Belyi function), and the embedding of T in C* can be seen as the 
preimage of M U {oo} by 7 tt, the three classes of vertices being the preimages of 0, 1 and oo 
respectively. So, the question of computing the embedding turns into that of computing ir T . 

The covering 7Tt is meromorphic, and in the case of the sphere it means that it has to 
be a rational function; computing it then means computing the location of its zeros and 
its poles, or equivalently the coefficients of its numerator and denominator. It is always 
possible to choose the embedding (via a Mobius transformation) in such a way that all 
these numbers are algebraic, and what we are looking for is their specification as roots of 
explicit polynomials with integer coefficients. This has been done in many cases before, cf. 
for example [5, 6, 13, 15], tracing back to the original work of Klein on the icosahedron and 
its link with fifth-degree equations [10]; a recent work with a similar goal as ours, but very 
different methods, is [3]. 

The same construction can be performed from a triangulation of any orientable surface, 
and leads to a Riemann surface Mt of the same genus together with a ramified covering 
7 Tt ■ M t —» C*; this surface in turn can be seen as a (smooth irreducible projective) algebraic 
curve, and Belyi’s theorem (see below for a precise statement) ensures that this curve is 
defined over Q. A similar goal as before is then to identify this algebraic curve explicitly, by 
specifying the coefficients of its equation as roots of explicit integer polynomials, as well as 
those of the covering. 

In this paper, besides the sphere we will only discuss the simpler case of genus 1, where the 
Riemann surface built from T can always be uniformized into a torus T = C/(Z + rZ) for some 
modulus r (depending on T ) in the complex upper-half plane. The previous problem becomes 
the exact computation of r and of the covering 7 Tt : T —> C*. Lifting 7Tt to the universal 


DESSINS D’ENFANTS FOR ANALYSTS 


3 


cover of T, this means that we are looking for an elliptic function rather than for a rational 
one as above, but apart from that the situation remains very similar in principle, in that we 
are looking for a finite collection of algebraic numbers. The higher-genus situation should be 
amenable to a similar treatment as the one we describe here, but the implementation would 
have to be significantly more complex. 

The usual ways that the question has been addressed previously are mostly algebraic in 
nature: one can write the conditions that a function has to satisfy in order to be a Belyi 
function as a system of polynomial equations, and try to solve it exactly from the start. This 
works well for small examples, but the combinatorial complexity grows very quickly and even 
moderately large cases, especially in positive genus, are beyond the reach of computer algebra 
systems. A recent proposal by Bartholdi et. al. [3] for instance quotes a computing time of 
15 minutes for a triangulation of the sphere with 15 vertices. 

Here instead we base our approach on a numerical point of view, and the central step is 
an implementation of the Newton algorithm to get an arbitrary precision approximation 
of the algebraic numbers we are looking for, initialized with a configuration obtained using 
circle packings; this can then be combined with a lattice-reduction algorithm to identify the 
corresponding integer polynomials. For comparison, for the same example as used in [3], 
getting 100 digits of precision for all parameters (which is more than enough to obtain an 
exact solution) takes less than a second. 

Our initial motivation for this work was related to probabilistic conjectures on the confor¬ 
mal structure of large random planar maps; testing these conjectures numerically without 
computing the embedding exactly is difficult, because in this setup approximations come 
from many sources: from the randomness of the map, from the fact that it is of finite size, 
and also from the embedding approximation itself, so determining the embedding exactly 
is of practical importance. Computing the embedding algebraically being beyond reach of 
the previous methods for the numbers of vertices that are relevant in this setup, we were led 
to look for a more numerical approach which is the focus of this paper; we will apply it to 
random triangulations in a subsequent work. 

The remainder of this paper is organized as follows: section 1 contains the formal definitions 
of the objects that we are using, section 2 gives a statement of Belyi’s theorem and is there 
mostly to keep the argument self-contained, and sections 3 and 4 describe our approach in 
detail and provide proofs of convergence in the cases of genus 0 and 1, respectively. Section 5 
then gives a few examples of application. 

The source code implementing the algorithms, and used to produce all the examples and 
pictures in this paper, is publicly available at the following address: 

http://github.com/vbeffara/Simulations 
1. Definitions 

1.1. Maps and hypermaps. Our starting point will always be a graph drawn on a given 
topological surface, but seen as a combinatorial object, i.e. up to homeomorphisms of the 
surface. All the following definitions are classical, but we still give them for sake of self- 
containedness. Much more can be found for instance is the book of Lando and Zvonkin [11], 
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Definition. A (finite) graph is a pair G = ( V , E) where V is seen as the set of vertices of G 
and E C V 2 as the set of its edges. We will always consider undirected graphs, for which E is 
symmetric: (x,y) G E <=$■ (y,x) G E, and assume the absence of loops, i.e. of edges of the 
form (x, x ) finking a vertex to itself. 

Definition. An embedding of an undirected graph G = {V, E) in an oriented surface M is 
the data of a collection (x v ) ve v of pairwise distinct points of M, and a collection ( / y x ,y)(x,y)eE 
of continuous simple curves on M, with p xyj : [0,1] —» M satisfying 7(0) = x and 7(1) = y 
and the symmetry condition 7 x , y {t) = J y , x (1 — t), and such that curves corresponding to 
distinct edges have disjoint images except possibly for their endpoints. The embedding is 
proper if every connected component of the complement of the union of all the images of the 
7 x>y (in other words, each face ) is homeomorphic to a disk. 

Definition. A map is an equivalence class of proper embeddings of a graph G in a surface 
M, where two embeddings are identified if they are conjugated by a homeomorphism of 
M. When M is the Riemann sphere C*, we will speak about a planar map ; when M is a 
two-dimensional torus, about a toroidal map. If each of the faces has exactly 3 edges along 
its boundary, the map is called a triangulation of M. 



Figure 1. A planar map with 5 vertices, 6 edges and 3 faces (including the outer 
one in this representation). 

A map is in fact a combinatorial structure, in the sense that there are only finitely many 
maps consisting in embeddings of a given finite graph in a given surface; it can be specified 
by ordering the neighbors of each of the vertices of the graph in a cyclic order. This remark 
allows for the definition of the dual map of a map: its (dual) vertices are in bijection with 
the faces of the map; two dual vertices will be declared adjacent if the corresponding faces of 
the initial map share an edge, and the neighbor ordering in the dual map is then given by 
the cyclic order around the corresponding face of the primal map. 

A very related structure, which makes computer implementations easier to manage in 
practice and is of algebraic relevance, is the following: 

Definition. Given a positive integer n, a hypermap of size n is a triple ( a , a. ip) of permutations 
of a set of n elements satisfying the compatibility relation amp = id. 

This can be seen as a generalization of the notion of map: indeed, restricting to the case 
where the cycles of a are all transpositions, one can interpret the domain of the permutations 
as the set of all the half-edges of the map, the cycles of a as its vertices (reading the half-edges 
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incident to that vertex in counterclockwise order), those of a as its edges (which are pairs of 
half-edges) and those of ip as its faces (reading the half-edges emanating from the vertices 
along that face in direct order) — see Figure 2 for an example, which is much clearer than 
any formal description would be. 



Figure 2 . Numbering of the half-edges of the map from Figure 1 to obtain its 
hypermap representation: in this case, we obtain a = (02)(164)(3 5 10)(78)(9 11), 
a = (0 1)(23)(45)(67)(89)(1011) and p = (043)(1 2 1097)(568 11). Following the 
usual convention, the label of each half-edge is to its left. 

The maps given as examples in this paper (for instance in Appendix 5) are described in 
this form, with the n-element set chosen as {0,... ,n — 1}. The triple corresponding to the 
dual map is then (ip -1 , a -1 , a~ l ) — and obviously here a~ x = a, but this way of denoting it 
makes the notation clearer. Conversely, hypermaps can also be seen as a particular case of 
maps: 

Definition. A graph is called bipartite if its vertex set can be partitioned into two disjoint 
subsets, in such a way that each vertex in one of them is only adjacent to vertices in the 
other one. A map is bipartite if the associated graph is bipartite; if the map is proper, this is 
equivalent to saying that its faces all have an even number of edges along their boundaries. 

It is customary to refer to the partition of the vertices of a bipartite map into black and 
white vertices. Such a map can then be encoded as a hypermap, where the cycles of a (resp. 
a, (p) correspond to the black vertices (resp. white vertices, faces). In this case, the cycles of 
(p have a length equal to a half of the number of edges on the boundary of the corresponding 
faces of the map. If every white vertex has degree 2, then the black vertices themselves form 
a graph with the same hypermap representation. 

One can compose these two constructions, starting from a map, seeing it as a hypermap 
where the cycles of a have length 2, and then seeing this hypermap as a bipartite map where 
the black vertices correspond to the vertices of the initial map and the white ones (with 
degree 2) to its edges. More graphically, this is equivalent to adding one vertex on each edge 
of the initial map, as shown in Figure 3. 

The last identification that we will need is with particular triangulations: 

Definition. A triangulation of a surface is called tripartite if its vertex set can be partitioned 
into 3 disjoint subsets, in such a way that each of its faces has one vertex of each of these 
subsets along its boundary. This is equivalent to saying that its dual map is bipartite, and 
that all its vertices have even degree. 
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Figure 3. The bipartite map corresponding to the hypermap in Figure 2. 

The three vertex sets will be represented here as black, white and red; in [ 1], vertices 
of the third kind are represented as asterisks. To each tripartite triangulation correspond 
3 bipartite maps, obtained by keeping the vertices in two out of the three subsets in the 
definition and the edges connecting them. Conversely, every proper bipartite map can be 
obtained this way, and another way of stating this is that every proper bipartite map can be 
completed into a tripartite triangulation of the same surface by adding a vertex inside each 
face and connecting it to the original vertices on the boundary of that face. 

Combining the previous remarks, every (proper) planar map can be refined into a tripartite 
triangulation by adding one vertex on each of its edges and one vertex inside each of its 
faces, connecting them in the natural way. If (a ,«, cp) is the hypermap corresponding to the 
initial map, then the cycles of each of these permutations are in bijection with the vertices 
in one of the three subsets in the partition of the tripartite triangulation. We will refer to 
this triangulation as the tripartite refinement of the map — see Figure 4. In the diagrams 
of Section 5, the edges incident to red vertices are not represented in order to make the 
combinatorics more readable. 



Figure 4. The tripartite refinement of to the planar map from Figure 1. One can 
check that the degree of each red vertex is equal to twice the number of edges of the 
corresponding face in the initial map. 
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1.2. Constellations. Our main objects of study in this paper are holomorphic maps defined 
between Riemann surfaces, but it is convenient as a way of encoding Belyi maps to introduce 
a more discrete kind of object. In this whole section, M will denote a fixed Riemann surface. 
Our first definition is that of a decorated point of M : 

Definition. A star on M is a pair z = (z, d ) where z £ M (the location ) and d is a positive 
integer (the multiplicity). 

The idea is that a star will specify the behavior of a map at its location (typically d would 
be the degree of ramification at z). A polynomial can be given as a finite sequence of stars 
located at its roots, and a rational fraction as two such sequences (for the roots and the 
poles). To match later with Belyi’s theorem, we will need just a little more information: 

Definition. A constellation on M is a triple C = ( Z,P,0 ) where Z = {z,}i^ nz , P = 
{P;}i^n p an d O = {oj}i^ no are three finite (potentially empty) sets of stars, which we 
will refer to respectively as the zeros , poles and ones of C. A constellation is non-degenerate 
if the locations of all involved stars are pairwise distinct (which we will always implicitly 
assume except otherwise mentioned). 

If / : M —y C* is a holomorphic covering of the Riemann sphere, it comes naturally with 
a constellation Cf = (Zf,Pf,Of) where Zf lists the preimages of 0, Pf those of 00 and Of 
those of 1, and the multiplicities are the corresponding orders of ramification. Under the 
assumptions that / only ramifies over {0,1, 00 }, / is known as a Belyi function, and Cf 
characterizes / uniquely; this relation can be made much more explicit in a few cases, which 
we describe now. 


1 . 2 . 1 . Polynomials on C. If / £ C[X] is a polynomial with complex coefficients and degree d, 
it is specified uniquely by the collection of its roots and their multiplicities, and its leading 
coefficient: if its constellation is given as 

C f = ({(^,d i )},{(oo,d)},{(o i ,<)}) 

(with Yhdi = d) then / can be recovered as 

f(z) = \l[(z-z t f 

where A is chosen to make the value of / equal to 1 at the points of Of. In this setup, / is a 
Belyi function if and only if fP d) — d, in which case 

f(z) - 1 = and f(z) = XdH(z-z i ) di - 1 l[(z- o l )<- 1 . 

The polynomial / is then known as a Shabat polynomial. 


1 . 2 . 2 . Rational fractions on C. Let / now be a (ramified) covering of the Riemann sphere 
by itself: it has to be a rational fraction with complex coefficients, and can be written 
as f(z) = p(z)/q(z ) where p (resp. q ) is a polynomial of degree d± (resp. d 2 ). Again 
Zf = {( Zi , di)} and Pf = {(pj. d) )} can be written explicitly in terms of the roots of p and q 
respectively, and 00 will occur in Zf, in Pf or in neither according to the relative values of d\ 
and d 2 . f can then be similarly reconstructed as 


/0) = A 


n (z - Ziff 

II(* - PiY x 


(1) 
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(where the products omit the potential term at oo, and where as before A is chosen so that / 
takes value 1 at the o t ). Once more, the constellation Cf corresponds to a Belyi function if 
and only if the sums of multiplicities are the same for Zf, Pf and Of. 


1.2.3. Functions on a complex torus. The case of genus 1 can also be made quite explicit. 
We will assume in this whole section that M is the complex torus with periods 1 and r where 
Q(t) > 0, in other words 

M = C/(Z + rZ). 

We will abuse notation by writing the locations of stars as complex numbers. Here, M comes 
with a natural covering 7r : C —>• M and if / : M —>• C* is holomorphic, it can be lifted as a 
periodic meromorphic function / : C —> C* (an elliptic function). 

The situation is a little more rigid than before, in the sense that some relations are 
automatically satisfied by the zeros and poles of /. More specifically, if Zf = {( Zi,di )} and 
Pf = {( Pi,<?i )} then 

(2) ^2 d i = 22 d 'i and 22 diZi = 22 d i Pi 


(the second relation being meant mod Z + rZ). The reconstruction of / from these can be 
made in a similar fashion as in the case of rational fractions: if the pi and are chosen 
in such a way that Y!= Y d'iPi (this time as complex numbers in C), then / can be 
obtained as 


( 3 ) 


f ( z ) = A 


nc 

nc (z-Pi) c 


in terms of the Weierstrass function ( with (quasi-)periods 1 and r. It is also the case 
that / can be written as a polynomial in the Weierstrass function p with the same periods 
and its derivative p' , or equivalently as a rational function on the associated elliptic curve, 
but although that is the “right” way to reconstruct / in terms of algebraic geometry, the 
coefficients of this polynomial depend on the constellation in a way that is less transparent 
and does not exhibit the same similarity with the previous two cases. 


1.2.4. A few additional remarks. In the case of Belyi functions, it is always the case that the 
sums of multiplicities is the same for Zf, Pf and Of, being equal to the degree of the covering; 
and the genus of the surface M is related to these multiplicities by the Riemann-Hurwitz 
formula. It is therefore natural to define the following: 

Definition. A constellation C = (Z, P, O) is balanced if it satisfies 

E d = E d = E d=:N(C). 

(w,d)£Z (w,d)eP (w,d)eO 

If C is balanced, its genus is given by 

g (C):=l + NiC )-* (ZUPU °) . 

2 

With these notations, if / : M — >■ C* is a Belyi function, then its constellation Cf is 
balanced and its genus it that of M. As we saw previously, if this genus is 0 or 1 , it is 
then possible to explicitly reconstruct / from C. Besides, if (w,d) G Of, then d gives the 
degree of ramification of / at w. This can be used in the other direction to characterize the 
constellations corresponding to a Belyi function: 
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Definition. Let C = (Z. P , O) be a balanced constellation of genus 0 or 1: we will say that it 
is exact if there exists A G C satisfying the following conditions. Let / be defined according 
to either (1) or (3) (depending on the genus — and assume that (2) is satisfied in the case of 
genus 1). Then for every (w, d ) G O : 

• f{w) = 1; 

• for every k G {1,. .., d — 1}, f( k \w) = 0. 

Such a value of A is obviously unique, and will be called the canonical normalization of C. 

1.3. Drawing a constellation. Let M be of genus 0 or 1 and C = (Z, P , O ) be an exact 
constellation on M; let / : M ^ C* be the corresponding covering. This setup allows to 
draw two natural, related structures on M. 



Figure 5. Uniformization of a periodic triangulation of 90 vertices. The color 
corresponds to the sign of the imaginary part of the covering map. 

1.3.1. As a triangulation. First, one can look at the tripartite map on M with vertices at 
ZUPUO, and edges given as the preimages of the three intervals (—oo, 0), (0,1) and (1, +oo) 
of R by /. This is a triangulation of M, and its dual is a bipartite map with all vertices of 
degree 3. Most of the pictures in this article are drawn in this setup (see e.g. the first page 
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for an example in genus 0, or Figure 5 in genus 1), and the faces are colored according to the 
sign of the imaginary part of / on them. 

1.3.2. As a hypermap. Another way to draw C is as a hypermap, where Z consists in the 
black vertices, P in the white vertices, and the edges are the preimages of the segment [0,1] 
by /. The combinatorics of the hypermap (the triple of permutations (<r, a, <f) in the language 
of the previous section) is equivalent to that of the triangulation above, and is completely 
determined by the constellation. It is this particular way of drawing this hypermap on M 
which is usually referred to as a dessin d’enfant. Even though the information is theoretically 
the same as before, it is computationally more difficult to obtain visually satisfactory pictures 
this way. 


2. Belyi’s theorem 

We give in this section an extremely limited introduction to Belyi’s theorem, for the benefit 
of readers who are not familiar with the topic. For much more, we refer to [9] and references 
therein. 

2.1. Preliminary remarks. For fixed balanced multiplicities, the space of constellations is 
of complex dimension #(ZUPU0). The number of equations to be satisfied for a constellation 
C to be exact is counted a bit differently depending on the genus: 

• If g — 0, it is the number of conditions coming from the definition, which is equal to 
N(C) - 1; 

• If g — 1, the relation (2) needs to be counted as well and the total number is IV (C). 

The difference 5 between the number of variables and the number of conditions is then equal 
to 3 in genus 0, and to 0 in genus 1; in particular it depends only on the genus. 

Before we can make use of this computation, one remark is in order. In the general setting, 
automorphisms of M act in a natural way on constellations on M by mapping locations while 
preserving multiplicities. This preserves the property of being exact; it is natural to identify 
two constellations conjugated by such an automorphism, or at least to classify them up to 
automorphism. 

In genus 0, the group of automorphisms of the Riemann sphere has complex dimension 3, 
which is equal to 6. This means that we should expect the number of exact constellations of 
given multiplicities, up to automorphisms of M, to be finite. 

In genus 1, the group of automorphisms of M has complex dimension 1 (in the generic case 
it consists purely of translations), which is one more than 5. This means that in the general 
case we should expect the existence of no exact constellation with given multiplicities, but 
that for finitely many values of r there should exist finitely many exact constellations. 

Of course, all the preceding remarks are at the heuristic level and making them formal 
would imply controlling degeneracies and genericity, which would likely be quite difficult to 
do at this point. Nevertheless, as what follows will show, they do give the right predictions. 

2.2. Analytic statement. Let T be a map with triangular faces, of genus 0 or 1 . It can be 
made into a Riemann surface Mt by gluing together equilateral triangles according to its 
combinatorics; and this surface can be uniformized to either the Riemann sphere, or to a 
complex torus with periods 1 and tt, where ly is uniquely determined (up to ,SX'_>(Z) action) 
by T. It is a natural question, and the origin of the work presented here, to in the latter 
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case determine the value of tt from the discrete data of the map T, either exactly or at least 
numerically (which, if done with enough precision, suffices to obtain an exact value). 

If the vertices of T all have even degree (which is in particular the case if T is associated to 
a hypermap), they can be seen as stars with multiplicity half of their degree and it is always 
possible to split them into three disjoint sets (or types) Z, P and O in such a way that each 
of the faces contains one vertex in each of these sets. One can then use the uniformization 
obtained above to map Mt onto the Riemann sphere, each triangle being mapped to either 
the upper or lower hemisphere according to the order in which the vertex types occur on its 
boundary. This construction leads to a covering f T : M T —» C* which is a Belyi function, and 
( Z , P , O) is then the exact constellation associated to fr ; see Figure 5 for an illustration in 
the case of genus 1. 

This gives a justification to the predictions at the end of the previous section: since each 
exact constellation gives rise to a triangulation, it can be obtained (up to automorphism) 
from that same triangulation. In turn, since the number of triangulations of a given size is 
finite, this implies that the number of exact constellations of given multiplicities, counted 
up to automorphism, is itself finite, and so is the set of values of r such that C/(Z + rZ) 
supports an exact constellation of a given size. 

This approach is very satisfactory at the theoretical level; unfortunately, it doesn’t lend 
itself very well to explicit computations, and performing this program numerically in a direct 
way would be difficult. 

2.3. Number-theoretic statement. There is a deep connection between the complex 
structure introduced above and the arithmetic properties of the underlying Riemann surface, 
having its origin in the following equivalence: 

Theorem 1 (Belyi). Let C be a non-singular algebraic curve: then there exists a ramified 
covering C —>• C*, ramified only above {0,1, oo}, if and only if C is defined over Q. 

We will mostly be concerned with the “only if” part in what follows; notice that when 
we build a Riemann surface from gluing equilateral triangles, it automatically comes with a 
constellation (Z, P, O) and such a covering, namely the function / constructed above, sending 
the points in Z (resp. P, O ) to 0 (resp. oo, 1) and unramified above the other points. 

What it means for us is that the torus C/(Z + rZ) built from a triangulation of genus 1, 
seen as an elliptic curve, is defined over Q. In particular, the modulus r can theoretically be 
computed explicitly. In practice, this seems to be possible only in very few cases, and in each 
of them, requires the use of particular symmetries of the triangulation. 

3. Numerical computation of Belyi maps: genus 0 

We now turn to the main point of this paper, namely a semi-numerical strategy to exactly 
compute the branched covering associated to a given hypermap. As argued above, this can 
be reduced to determining the exact constellation associated to the map, which lends itself 
to numerical approaches. We start with the case of hypermaps drawn on the sphere, which 
benefits from lighter notation while retaining most of the features of our approach. 

Fix three tuples of positive integers, d° = (d°)i^ no , d°° = (d?°)i<^noo and d 1 = (d})i^ ni 
with the same sum N , satisfying the genus-0 condition 


no + Uqo + Tii — N + 2 
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and let £ = <Z ( p,d° o ,d 1 be the space of all non-degenerate constellations on C* with signature 
(d°, d°°, d 1 ) (which are all balanced and of genus 0). £ is a complex manifold of dimension 
N + 2. If C = {Z, P, O ) G £, let 

n (. z-w) d JJ (z-w)~ d 

(w,d)€Z (w,d)£P 

be defined as above (omitting the star at oo in the products, if any), and consider the map 
$ : C x £ —» C N defined by 

$(A,C) := (/(oi) - l,/ / (oi),...,/ (d i“ 1) (o 1 ),/(o 2 ) - l,/ / (o 2 ),...,/ (d "i“ 1) (o ni )), 

where to lighten notation we let / = f\ t c and C = ( Z,P,0 ) with O = ((oj,d|)). The 
constellation C is exact if and only if there exists A G C such that 4>(A, C) = (0,..., 0). 

To normalize the map embeddings, let €* be the set of all constellations C = (Z. P. O) G £ 
such that z\ = 0, pi = oo and Oi = 1. (£* is a manifold of complex dimension N — 1 which 
can be seen as an open subset of C iV_l by listing the locations of all the stars except for z i, 
Pi and o 1 . The restriction <f>* : C x £* —> C N can thus be seen as a map from an open subset 
of C N to C N , which is clearly analytic in all its variables. 

From the previous discussion, the set of preimages of (0,..., 0) by <b* is finite, and each 
of its elements corresponds to an exact constellation and to the associated hypermap and 
covering. From this, the general structure of the algorithm is rather clear: starting from a 
tripartite hypermap M and the associated triangulation of genus 0, to find its constellation 
Cm we will perform the following. 

• Step 1: find an approximation Co of Cm', 

• Step 2: refine the approximation via an iterative scheme, using Co as a starting point; 

• Step 3: verify that the solution obtained is the right one (topologically); otherwise go 
back to step 1 with better precision; 

• Step 4: identify the coefficients as algebraic numbers; 

• Step 5: verify that the solution obtained is the right one (algebraically). 

We now describe each of these steps in some detail, focusing more on the first two which 
constitute the main contribution of the section. 

3.1. Step 1: Approaching the complex structure. To initialize the algorithm, we need 
to get an approximation of Cm- This can be done by approximating the complex structure 
derived from M, and there are several options available here. From the complex analytical 
point of view, one natural possibility would be to start from an arbitrary embedding of M 
into the sphere with straight edges, and to identify the uniformizing map as a solution to 
a Beltrami equation with a Beltrami differential taken to be constant on each of the faces. 
This has been implemented in a few cases but is quite involved [8], so we chose a simpler way, 
following [7]. 

Given a triangulation T of the sphere, there exists a circle packing with the combinatorics of 
T, namely a collection (C v ) of circles indexed by the vertex set of T, such that the associated 
disks have disjoint interiors and such that C v and C w are tangent if and only if v and w are 
adjacent vertices of T. The circle packing is unique up to Mobius transformations, and it 
provides an embedding of T in which the vertices are mapped to the centers of the circles 
and the edges are each the union of two circle radii; in particular, if T is tripartite it provides 
also a balanced constellation C cp (T). For general reference on circle packings, we refer the 
reader to [14] 
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Our first approximation of Cm will simply be C cp (M). In practice, it turns out to often be 
good enough for our purposes, but it is not always the case. To get a better one, we follow 
the strategy of [ ] and refine M in the following way. Given a triangulation, one can replace 
each of its faces by 4 new triangles by inserting a vertex at the mid-point of each of its edges 
— see Figure 6. Iterating the procedure starting from M 0 := M, one gets a sequence (M n ) n ^ 0 
of triangulation, which are all tripartite. The main result is then the following: 





Figure 6. Successive refinements of a triangulation. 

Theorem 2 (Theorem 4.8 in [ ]). Under appropriate normalization, C cp (M n ) converges to 
Cm as n ^ oo. 

This is exactly what we needed: we now have a sequence of constellations which approaches 
Cm and will serve as starting points for the iterative method in the next step. We first run steps 
2 and 3 starting from C cp (M 0 ) = C cp (M ), then if that fails from C cp (Mi), C cp (M 2 ), ... until 
the first one that succeeds — which is guaranteed to happen after finitely many subdivisions 
by Theorems 2 and 3. 

3.2. Step 2: Newton’s method. This part is quite standard, but as far as we can tell it 
has not been implemented in this setup before: one can apply a multi-dimensional Newton’s 
method to build a sequence of constellations which converges to Cm, provided one starts close 
enough to it. 

More specifically, still identifying (£* with an open subset of CW -1 , let T : C x €* —* C x (£* 
be defined by 

*((A,C)) := (A,C) - J t ((A,C))-'4((A,C)) 

whenever the Jacobian is nonsingular at (A. C). The main statement of this section is the 
convergence of Newton’s method: 

Theorem 3. Let A m be the canonical normalization of Cm-' there exists a neighborhood U of 
(A m,Cm) in C x (£* satisfying T(N) C U, on which is everywhere nonsingular, and such 
that, whenever (A ,C) G U, the sequence defined inductively by Xq = (A ,C) and x n+ i = T(x n ) 
converges to (A m,Cm) as n — »■ oo. 
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It will be more convenient to set up the computations in the proof in terms of logarithmic 
derivatives, so let g\,c(z) be a determination of log(/(A,c)(^)) that is jointly continuous in 
(A,C) in a neighborhood V of (A m,Cm) and in 2 in neighborhoods of each of the o iy and 
chosen such that g(\ M ,c M )(°i ) = 0 for ever Y * £ {1,..., ni}. Let A be the counterpart of d> 
defined in V in terms of g instead of /, namely: 

A(A,C) := (g(oi),g'(oi ),..., ^ _1 ) (oi), g(o 2 ), g'(o 2 ),..., g {d ^~ l) {o ni )). 

Lemma 4. There is a constant C £ C \ {0}, depending only on the signature (d°, d°°, d 1 ), 
such that the Jacobian determinant of A at (A m,Cm) Is given by 




2^i<j^no 


( z j z i) U 


2^i<j^n c 




A 




\d} 


where d] := d\ — 1 + Si- In particular, det Ja((\m,Cm)) 7 ^ 0 . 

Proof. We first need to compute the entries of the Jacobian matrix Ja((Am,Cm))- The 
iterated derivatives of g are given for k > 0 by 


(-l ) fc+1 

Jk^Ty. 


no 

9 W M = £ 
1=1 



E 


df 


■rf ( w ~ z i) 


.\k 


so the partial derivatives are all explicit: 


dxg^M = i<5», 


9» J 9 < ‘ l (o i )=S ( ‘ +1| (o i )^, 


d Zj g ik \oi ) 


dfl(-l) k+l k\ 


Zj ) k + l ’ 


dfjg {k) (°i ) 


d°°(-l) k+1 k\ 

(oi - /i ) fc+1 ’ 


In addition, at the point (A m,Cm) we know that g^ k \oi) = 0 as soon as k < d], so up to 
terms depending only on the signature (and which are products of (— l) k , dj -s and factorials), 
the Jacobian determinant is proportional to 


l/A 

(01 -Z2)- 1 

o--- 0 

(01 - z 2 )- d ' 

1/A 0-- 

(o 2 - Z2Y 1 

0 

(02 - Z 2 )~ d 2 • • 

1/A 0-- 

(On% ~ Z 2)- 1 

0 

Om - z 2 y d1 ^ 

(Ol - ZnoY 1 

(01 - A ) -1 

■■■ (01 - z„ 0 ) _di 

(Ol- A)”^ 

(o 2 ~ z no ) 1 ••• 

(o 2 - A ) -1 

(02 - Z n 0)-^ • • 

(02 - A) _d = • 

(°ni — z n 0 ) 

* {0 ni ~ / 2) _1 

(°n 1 — z n 0 ) ni 

(■ o ni ~ / 2 ) -41 

(Ot-fnj - 1 

0 

••• Dl-inco)”^ 

0 

(02 ~fn J _1 ■■■ 

0 

(02 -AJ -4 •• 

1 

(°m — fn 00) 

0 

(°nr - /nj "* 1 

0 

0 

0 

0 

0 

0 

1 


where the columns are the components of A and the lines correspond to partial derivatives with 
respect, in order, to A, the (zf )^ 2 , the (ff )^2 and the (oj)^ 2 (remember that we normalized 
all constellations to have z\ = 0, Oi = 1 and f\ = 00 , so they do not appear as variables here, 
and that this makes the matrix square). Developing the determinant along its last (711 — 1) 
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lines shows that it is equal to 


l/A 0••• 

0 

l/A 0--' 

0 

l/A 0• • ■ 

0 

'cP 

1 

O 

to 

1 

(01 - z 2 ) _ 4 

O 2 - 22 ) _1 

( 02 _ 22 )i-4 ... 

On, “ 220 1 

On, - 2 2 ) 1-4 ‘ 

Ol Z no ) 

Ol - ^no) _ 4 

O 2 — ^no) 

02-2n„) 1 -4 ••• 

(°n i z no) 

(On, - Z n0 ) 1_d h 

(Ol - / 2)- 1 • • • 

Oi - /2)“4 

O 2 “ T 2 ) 1 

(0 2 - /2) 1 - 4 • • • 

Om - /fi” 1 

(On, - f2f~< 

(01 - /noo) _1 • • • 

Ol - fn 00 ) “4 

(02-/noJ _1 ••• 

(02-/noo) 1_ 4 ••• 

(°ni ~ /noJ -1 ' ' ' 

(On, - fnj 1 ^ 

which can now 

be identified 

as a rational fraction in the z*, 

fi and Oi. 



The denominator is easier to compute first: it just consists in the product of the terms 
(oj — z 3 ) k and (o; — f :] ) k with the highest power appearing in the determinant, namely 
dj — 1 + 5} = d] (which is the width of the matrix block corresponding to of). Multiplying 
each line by those factors which appear in it, and the first one by A, one gets a matrix A with 
polynomial entries in the Zi, fi and oy let P be the determinant of that matrix. To compute 
it, we now forget about the origin of the matrix and note that if any two among {z^} U {/,} 
are equal, then two of the lines of A are equal, hence P is divisible by the product of all the 
differences between these values. 

If on the other hand o, t = Oj for some i j, then two blocks of vertical lines become equal, 
so P is divisible by some power of (o l — Oj). Getting the value of the power is a little bit 
more involved but again can be done very explicitly. Start with the following identity (where 
a ^ x and k G Z + ): 

. 11a: x k ~ k x k 

(4) -= - + — + ••• H j— + —: -r. 

a — x a a z a K a K (a — x) 

This can be applied to get an expansion of the first column in the block corresponding to of 
for the second line, letting a = (oj — z 2 ), x = (o* — Oj ) and k = dj — 1 + Sj gives an expansion 
of (Oj — z 2 )~ k in terms of the entries in the same line in the block of Oj, with coefficients 
depending only on (o, — Oj ) and hence being the same across all lines, plus a remainder where 
(oj — Oj) k is in factor: more explicitly, 

1 1 \ Qi ~ 0j \ | | 

Oj — z 2 Oi- z 2 (Oi - z 2 ) 2 {Oi - z 2 ) k ( Oi - z 2 ) k (oj - z 2 )' 

The next columns can be expanded in a similar way, where the expansion of the Gth 
column in the block of o 3 uses the last (dj — 1 + Sj) — (£ — l) column of the block of o t and 
the remainders in the first (£ — 1) columns of the block of Oj, with coefficients depending only 
on (oj — Oj) and a remainder term where the same power of (o l — Oj) as before factors out. 
These expansions are obtained from successive derivatives of (4) with respect to x; the first 
one being for instance 

1 1 2x (k — l)x k ~ 2 kx k_1 x k 

(a — x) 2 a 2 a 3 a k a k (a — x) a k (a — x) 2 

Overall, the power of (oj — Oj) in P obtained this way is exactly the product of the widths of 
the corresponding blocks in the matrix. 

To summarize, P is divisible by the numerator in the statement of the theorem. Matching 
the degrees shows that the ratio between det Ja and the formula in the statement of the 
lemma is a constant, which is the main claim. The fact that the determinant is not zero is 
then a direct consequence of the fact that Cm is non-degenerate. □ 
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Proof of Theorem 3. First, we can replace / with e 9 in the expression for <3>, and compute 
the partial derivatives appearing in J<\, in terms of those appearing in L\. By repeated use of 
the chain rule, commutativity of derivatives and the fact that Cm is exact, we get for example 


d Zi [ffx,c) {j \ok)] = (. f\,cd Zi g\,c) U) (ok ) = ^2 = d Zi [(g\,c) U \o k )] 


e=o 


at the point (A m,Cm) as soon as j < d\. Similar computations in the other variables show 
that in fact J\ and are equal at the point (A m,Cm) and in particular, J$((Am,Fm)) is 
non-singular. 

The remainder of the proof is a completely standard application of the usual Newton 
algorithm: it suffices to use the smoothness of <I> in all variables to show that J$((A, $)) is 
non-singular in a neighborhood of (A m,Cm), and to expand T at (A m,Cm) to obtain a bound 
of the form 

ll'^(('V£)) — (Am,Cm)|| C1|(A,C) — (Am,Cm)|| 2 

for (A ,C) close enough to (A m,Cm)- □ 


To summarize our construction so far, we are essentially done at the theoretical level: there 
is a neighborhood U of our point of interest (Am, Cm) from which Newton’s algorithm is guar¬ 
anteed to converge, and iterating the subdivision in the first step will bring the approximation 
within U after finitely many steps. In addition, since the convergence is quadratic, one can 
iterate a computer implementation until two successive values are indistinguishable within 
machine precision (or chosen extended precision) at very little cost, meaning that we can 
get a numerical approximation of (Am, Cm) with arbitrarily chosen precision in a reasonable 
computing time. 


A significant issue in practice though is that none of the two steps is quantitative. The speed 
of convergence of circle packing embeddings to the uniformizing map is not well understood 
(although our simulations as well as the numerical experiments at the end of [ ] suggest that 
it should be polynomial in the diameter of the smallest circle and exponential in the number 
of refinements), and the neighborhood U we would be able to explicitly construct by keeping 
track of all constants implicit in the proof above would certainly be much smaller than the 
basin of attraction of (Am,Cm)- 

What this means is that the number of subdivisions we should perform in step 1 to be 
certain to have convergence would make that step computationally unfeasible. We now turn 
to a way around this issue. 

3.3. Step 3: Topological verification. As it turns out, in practice the domain of attraction 
of Newton’s method in the case we are interested in seems to be quite large, at least much 
larger than continuity arguments for <f> would predict. A natural procedure is therefore to 
start Newton’s algorithm from successive iterations of step 1, starting in fact with no iteration 
at all, and in each case to see if the iteration converges or not (quadratic convergence when 
the iteration is successful means that recognizing convergence is very quick). 

If the iteration does not converge, or converges to a singular constellation, then we simply 
refine once more. If it does converge to a non-singular pair (A^, C^), then we found a 
constellation with the right signature. It remains to check whether it actually is the one we 
were looking for (and subdivide once more if it is not). 
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One can simply do it visually, from a picture of the sign of the imaginary part of the 
function fx^fioo like those in the figures of this paper, recovering the triangulation and just 
checking that it is isomorphic to M. The same thing can probably be automated, although it 
is not at all clear which algorithms could be proved to work; constructing a region adjacency 
graph from the picture (thus recovering the dual graph of the triangulation) would be an 
option, but we did not attempt to implement it. 

In practice, except on specially tuned cases designed to test the numerical stability of the 
whole method (similar for instance to the map shown in Figure 7 with vertices of very large 
degree), the limit is almost always the right one — though of course this statement is not of 
a mathematical nature and is merely an empirical observation. 

3.4. Step 4: Lattice reduction. At this point, we obtained a numerical approximation of 
(Am, Cm) with the normalization chosen in such a way that the function f\ M ,c M h as a zero f d 0, 
a pole at oo and takes the value one at 1. We know in advance that with such normalization, 
the locations of the other zeros, poles and ones are all algebraic numbers: indeed, there is a 
choice of normalization such that all locations are algebraic and mapping it to our preferred 
one can in turn be done by applying a Mobius map with algebraic coefficients. 

This means that for each of those algebraic numbers, we are able to obtain an approximation 
to any precision that we want, in a reasonable computing time. Identifying the minimal 
polynomial of an algebraic number given such an approximation is a much-studied question, 
and there exist a variety of classical integer relation algorithms to do it, based on lattice- 
reduction methods. Going into a detailed description of such algorithms is besides the point 
of this paper; several implementations are freely available, and the results presented below 
were obtained using one of them (specifically FPLLL, see [2]). 

While the empirical observation is that these methods work, again a comment of a more 
theoretical nature is in order: even though the numbers that we are interested in are all 
algebraic, the bounds on their degree that one can derive from the proof of Belyi’s theorem are 
enormous and the precision approximation that is required for lattice reduction to provably 
find the right solution is therefore enormous as well, to the point that making those bounds 
quantitative is likely to be of little practical use. 

A more interesting question is that of the choice of normalization. The one we chose was 
convenient from the implementation point of view, as it made the formulas in Step 2 explicit, 
but there is no reason why it would lead to the algebraic numbers of the lowest possible degree, 
thus compounding the previous remark. Other choices are of course available: for instance 
one might want the sum with multiplicity of all zeros to be 0 (to get a vanishing coefficient 
in the numerator of f\ M ,c M )- More convincingly, if M has a non-trivial automorphism group, 
so does its constellation, and one may want (some of) the corresponding automorphisms to 
be affine transformations of the plane. 

This however has little impact on the implementation, as one can always go from one 
normalization to another after having obtained the approximation in Step 3. 

3.5. Step 5: Algebraic verification. The last validation step is of a purely algebraic nature: 
given the list of locations as roots of integer polynomials, verify that the constellation they 
form is indeed exact, in which case we have achieved our programme of computing the Belyi 
function associated to the map M explicitly. This can be rewritten as a collection of algebraic 
equations that they must satisfy, and can be checked explicitly either by hand (for smaller 
cases) or by a computer algebra system; again, if the validation failed, it means that the 
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Figure 7. Uniformization of a triangulation of genus 1 with one vertex of high 
degree (which can be used as a test-case for numerical stability of the algorithm). 

polynomials obtained in the previous step are erroneous, and one can re-run the lattice 
reduction from higher precision approximations, with the guarantee that after finitely many 
round trips the right solution will be found. 

Remark 1. This last step is usually presented as the starting point of the computation of 
Belyi functions: namely, starting from a hypermap, obtain a system of polynomial equations 
in the location of the ramification points, and then use elimination theory to solve the system 
(usually using Grobner bases). This works well in practice, but does not seem to extend well 
to higher genus beyond very small maps. 

4. Numerical computation of Belyi maps: genus 1 

We now turn to the case of triangulations of the torus, and to elliptic Belyi functions. As 
mentioned earlier, the main structure of the construction is extremely similar to the one 
described in the previous section, and we will focus on the differences rather than giving a 
complete description, using the same notation for corresponding but slightly different objects 
where it doesn’t lead to confusion. 
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We will work on the space of all non-degenerate constellations of genus 1 and of a given 
signature defined on a complex torus; however, there is an invariant (the modulus of the 
torus) which we do not know in advance. This means that the objects we will work with are 
in fact triples of the form (r, A, C) where r is the modulus of a torus, A is a normalizing factor 
and C is a constellation in the torus T t = C/(Z + rZ), and having such a representation of 
the Riemann surface on which we are working will be very convenient. 

So again, fix three tuples of positive integers, d° = d°° = (d°°) Ki^noo and 

d 1 = (dj )with the same sum N, satisfying the genus-1 condition 


n 0 + noo + ni = N , 


and let X = r £ d o )d oo d \ be the space of all triples (r, A, C) where r 6 i, A £ C \ {0} and C is a 
non-degenerate constellations on T t with signature (ci°, d°°, d 1 ). X is a complex manifold of 
dimension N + 2. We will always see T as the quotient of the complex plane by the lattice 
Z + rZ, represent the locations of the stars of C as complex numbers and identify X as an 
open subset of CW" 1 " 2 ; if C = ( Z , P, O ) G £, let 


fr,\,C '■ Z A 


n cl - Pi) ir 


be defined as above, and consider the map $ : X —» C A ’ defined by 

$((tA,C)) := (/(oi) - l,/ / (oi),...,/ (d i“ 1 ) (oi),/(o 2 ) - l,/ , (o 2 ),...,/ (,i "i _ 1 ) (on 1 )), 


where to lighten notation we let / = / T; a,c- 

To normalize the map embeddings and ensure that / is doubly periodic with periods 1 
and r (or in other words, to have it defined on the torus T t ), let X* be the set of all triples 
(r, A, C) G X such that 

( 6 ) (J °i z * = ^2 d i° Pi = °- 

X* is a manifold of complex dimension N which can be seen as an open subset of C A "' by 
listing the locations of all the stars except for zi and p\\ the restriction : X* —> C iY can 
thus be seen as a map from an open subset of C iY to C ,Y , which is clearly analytic in all its 
variables. We are interested in finding one specific preimage of (0,..., 0) by T. out of finitely 
many. 


We are now exactly in the same situation as before, the only difference being the definition 
of /, so we will briefly review the relevant changes needed to the previous steps. 

4.1. Step 1: Approaching the complex structure. Here the construction starts from 
a triangulation T of genus 1, and again there exists a unique Ty £ I and an essentially 
unique circle packing on the torus of modulus Tt having the combinatorics of T. Equivalently, 
the universal cover T of T is a doubly periodic triangulation of the plane, and there is an 
essentially unique locally finite circle packing of the plane with the combinatorics of T ; this 
circle packing is automatically doubly periodic, and can be normalized in such a way that it 
has the two periods {1, tt }. 

We want to define the starting point for the Newton iterative scheme using tt as the 
modulus and the locations of the centers of the circles as those of the stars; and here as well, 
under successive refinements, these converge to the constellation corresponding to M. One 
slight problem is that there is no reason why the normalization condition (6) can be made 
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to hold by a suitable choice of embedding of a circle packing — but we can at least ensure 
that ( ^i z t = 0, which makes the embedding unique, and simply ignore the sum of the d^p, 
(which will tend to 0 as the triangulation is refined more and more). 

4.2. Step 2: Newton’s method. This part works exactly the same way as before, with 

T((r, A, O) := (r, A, C) - J*((r, A, C))- 1 $((r, A ,C)), 

and the only point to check is the non-singularity of the Jacobian matrix at the point 
(t m , A a/. Cm). This can be done in a very similar way as on the sphere, as far as the variables 
besides r are concerned: for fixed r the Jacobian determinant is an elliptic function in all its 
other variables, and it factorizes for the same reason as in the rational case. Alternatively, 
one can also get the non-singularity of the Jacobian by moving all the locations of the points 
closer and closer to the origin (while satisfying the constraints listed above), and noticing 
that the asymptotic behavior of the determinant, once properly normalized, is given by the 
case of genus 0. 

The variable r is a bit more problematic, because the derivatives in r of the functions 
involved in <f> are not as explicit as the derivatives in z. Instead, one way to go around the 
issue is to argue that a different way of proceeding, from the data of a toroidal triangulation, 
is as follows: take a large square ofiVxiV periods in the universal cover of M, and stitch it 
with a copy of itself along the boundary of the square to obtain a triangulation of the sphere. 
This triangulation can be uniformized (for instance using the method in the previous section) 
and the local behavior of the uniformizing function near the center of the square, as N —> oo, 
becomes periodic and asymptotic to the uniformizing map of M itself. In particular, the 
fact that the Jacobian is non-degenerate in genus 1 can be extracted from the corresponding 
statement in genus 0. 

Remark 2. While that last remark gives an alternative strategy to solve our initial problem as 
well, it has two main drawbacks: first, it is difficult to estimate the level of precision that one 
would be obtained as a function of N (presumably it would be polynomial in N ); second, and 
very related, is that to get the kind of approximation needed to apply the next step, the value 
of N would have to be taken so large that there would be no hope of actually implementing 
the programme. It would still be a possibility to obtain the starting point of step 2 in this 
way, thus replacing step 1; but we did not try this route. 

4.3. Steps 3, 4 and 5: Identification and validation. Here, not much needs to be 
changed at all except for notation; the questions raised, whether the successive objects that 
we construct correspond to the one we are looking for, are the same, and the method is 
exactly parallel to that used in genus 0. 

One point needs to be made though. In the sphere, the freedom in the choice of embedding 
meant that finding the constellation leading to algebraic numbers with minimal degree was an 
issue. Here, there is much more rigidity from the conditions (6), which itself is very natural, 
and there is no choice at all, beyond the usual discrete SL^ (Z) action, in the modulus r; 
indeed the elliptic curve defined by M is uniquely defined. This is quite convenient, especially 
since determining r was our primary goal in the case of genus 1. 

5. A FEW EXAMPLES 

The front page of this paper shows the Belyi function of the tripartite refinement of a 
uniformly sampled random triangulation of the sphere with 15 vertices. We validated the 
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algorithm of a few known cases, such as those listed in [1], but chose to focus on genus 1 for 
the examples of this section. 

We list here the Belyi functions obtained from all the triangulations of the torus with up 
to 3 vertices and all degrees at least equal to 3. In each case, the triangulation is refined into 
its tripartite refinement and we list the hypermap description of the original triangulation, 
a graphical representation of its combinatorics, and then the sign of the imaginary part of 
the Belyi function and in a few instances the exact value of the associated j -invariant (or its 
minimal polynomial in case that is more useful). For triangulations with 4 vertices we give 
one example in the same form, and only the graph of the covering map for the others. 


In several cases the original triangulation is itself tripartite, so the representation is not 
minimal; but the value of j is not affected by the refinement, and the dessin associated to 
the initial triangulation is a subset of that of the refined one. 

Each finite triangulation T has countably many coverings (Tj^) that are themselves 
finite triangulations, having as fundamental domains unions of finitely many copies of the 
fundamental domain of T arranged as a k x l rectangle. Those are omitted in the catalog 
below, and only the minimal one is listed. 


5.1. Size 1 . There is only one triangulation of the torus with one vertex, and its natural 
embedding is the usual triangular lattice composed of equilateral triangles: 



<x = (052143) 
a = (01) (2 3) (4 5) 
<p = (024) (135) 



3 = 0 


5.2. Size 2. In addition to the double coverings of the previous one, there are two minimal 
triangulations of the torus with two vertices: one is a refinement of the triangular lattice 
(with the same j invariant equal to 0) and the other is the face-centered square lattice. 
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a = (0785219104) (3116) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) 
(p = (026) (148) (3510) (7119) 





5.3. Size 3. There are 9 triangulations of the torus with 3 vertices and minimal degree at 
least 3 that are not coverings of previously displayed cases. They turn out to all have rational 
j-invariants. 
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a = (07171495126315104) (18 2) (1116 13) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (1213) (1415) (1617) 
if = (0 2 6) (14 9) (3 8 14) (5 10 13) (712 16) (1115 17) 



3 = 0 



a = (0 716 104) (19 12 63 1114 8 5 2) (13 15 17) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
f = (0 2 6) (14 8) (3 5 10) (712 17) (91413) (111615) 
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a = (0 714 104) (191715 13 8 5 2) (3111612 6) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
(^ = (02 6) (14 8) (35 10) (712 15) (91316) (111417) 



j = -3 072 



a = (0 5 13 10 6 2 1 714 8 4 3) (9 16 12) (111715) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
ip = (0 2 4) (13 6) (5 8 12) (710 15) (91417) (111316) 
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a = (0 7171510 4) (19161314 8 5 2) (311126) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
(^ = (02 6) (14 8) (35 10) (712 16) (91417) (1115 13) 




o = (0 51314 843) (171715106 2) (9161112) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
<p = (0 2 4) (13 6) (5 8 12) (710 16) (91417) (1115 13) 
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a = (0 716 8 2 11114 9 4) (3 15 12 6) (51713 10) 

a = (0 1 ) (2 3 ) (4 5 ) (6 7 ) (8 9 ) ( 1011 ) (12 13 ) ( 1415 ) ( 1617 ) 
(^ = (02 6 ) ( 1410 ) (3 814 ) (5 9 16 ) (712 17 ) ( 111315 ) 




(7 = (0 714 8 5 2191610 4) (3 1112 6) (13 1715) 
a = (0 1 ) (2 3 ) (4 5 ) (6 7 ) (8 9 ) ( 1011 ) (12 13 ) ( 1415 ) ( 1617 ) 
p = (0 2 6 ) (14 8 ) (3 5 10 ) (712 15 ) ( 91417 ) ( 111613 ) 
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a = (05131411843) (1715161062) (91712) 
a = (0 1) (2 3) (4 5) (6 7) (8 9) (1011) (12 13) (1415) (1617) 
(^ = (02 4) (136) (58 12) (710 14) (91116) (131715) 



1977326 743 



5.4. Size 4. With 4 vertices there are already too many triangulations to make listing them 
very useful. We still display one in detail, and will only show the picture for the others: 



This is the triangulation used as an example in [4, 12], where the elliptic curve is derived 
formally. It is defined over Q[v7] and one can check that the j-invariant obtained by our 
method is the right one (the discriminant of the polynomial above is 7 • (345 128) 2 ). Note 
though that the proof as detailed in [ ] is 4 pages long and moreover relies very strongly on 
the additional symmetries of the triangulation — in the picture above one can readily see 
that the embedding is symmetric under reflection by the line going through 0 and 1 + r and 
to another one orthogonal to it, which both correspond to automorphisms of order 2 of the 
triangulation. 

This is a strong indication that very small examples are computationally difficult to address, 
and slightly larger ones, or even cases of small size but no symmetry, are beyond these methods. 
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In comparison, the programme described here gets the exact constellation to machine precision 
(i.e., to within 10 -13 ) in less than a tenth of a second and producing enough digits to obtain 
the minimal polynomial for j takes of the order of 20 seconds on a standard laptop. 











WA 
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